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First order rotational perturbations of the Friedmann-Robertson- Walker metric are considered in 
the framework of the brane-world cosmological models. A rotation equation, relating the perturba- 
tions of the metric tensor to the angular velocity of the matter on the brane is derived under the 
assumption of slow rotation. The mathematical structure of the rotation equation imposes strong 
restrictions on the temporal and spatial dependence of the brane matter angular velocity. The study 
' of the integrable cases of the rotation equation leads to three distinct models, which are considered 

in detail. As a general result we find that, similarly to the general relativistic case, the rotational 
perturbations decay due to the expansion of the matter on the brane. One of the obtained consis- 
' tency conditions leads to a particular, purely inflationary brane-world cosmological model, with the 

, cosmological fluid obeying a non-linear barotropic equation of state. 

PACS numbers: 04.20.Jb, 04.65. +e, 98.80.-k 
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I. INTRODUCTION 



On an astronomical scale rotation is a basic property of cosmic objects. The rotation of planets, stars and galaxies 
inspired Gamow to suggest that the Universe is rotating and the angular momentum of stars and galaxies could be a 
result of the cosmic vorticity But even that observational evidences of cosmological rotation have been reported 
@-[D) they are still subject of controversy. 

From the analysis of microwave background anisotropy Collins and Hawking Q and Barrow, Juszkiewicz and 
Sonoda |Q have found some very tight limits of the cosmological vorticity, T i, s > 3 x 10 5 Th, where T b s is the actual 
rotation period of our Universe and Th = (1 ~ 2) x 10 10 years is the Hubble time. Therefore our present day Universe 
is rotating very slowly, if at all. 

From a theoretical point of view in 1949 Godel || gave his famous example of a rotating cosmological solution to 
the Einstein gravitational held equations. The Godel metric, describing a dust Universe with energy density p in the 
presence of a negative cosmological constant A is 



-(dt + e x dzf + dx 2 + dy 2 + ^e 2x dz 2 



(1) 



In this model the angular velocity of the cosmic rotation is given by uj 2 = Amp = —A. Godel also discussed the 
possibility of a cosmic explanation of the galactic rotation jg]. This rotating solution has attracted considerable 
interest because the corresponding Universes possess the property of closed time-like curves. 

The investigation of rotating and rotating-expanding Universes generated a large amount of literature in the field of 
general relativity, the combination of rotation with expansion in realistic cosmological models being one of the most 
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difficult tasks in cosmology (see M for a recent review of the expansion-rotation problem in general relativity). Hence 
rotating solutions of the gravitational field equations cannot be excluded a priori. But this raises the question of why 
the Universe rotates so slowly. This problem can also be naturally solved in the framework of the inflationary model. 
Ellis and Olive [0 and Gr0n and Soleng pointed out that if the Universe came into being as a mini-universe of 
Planck dimensions and went directly into an inflationary epoch driven by a scalar field with a flat potential, due to 
the non-rotation of the false vacuum and the exponential expansion during inflation the cosmic vorticity has decayed 
by a factor of about 1CU 145 . The most important diluting effect of the order of 10~ 116 is due to the relative density 
of the rotating fluid compared to the non-rotating decay products of the false vacuum . Inflationary cosmology 
also ruled out the possibility that the vorticity of galaxies and stars be of cosmic origin. 

The possibility of incorporating a slowly rotating Universe into the framework of Friedmann-Robertson- Walker 
(FRW) type metrics has been considered by Bayin and Cooperstock [|l2|, who obtained the restrictions imposed by 
the field equations on the matter angular velocity. They also shown that uniform rotation is incompatible with the dust 
filled (zero pressure) and with the radiation dominated Universe. Bayin has also shown that the field equations 
admit solutions for a special class of nonseparable rotation functions of the matter distribution. The investigation 
of the first order rotational perturbations of flat FRW type Universes proved to be useful in the study of string 
cosmological models with dilaton and axion fields fll4| . The form of the rotation equation imposes strong constraints 
on the form of the dilaton field potential U, restricting the allowed forms to two: the trivial case U = and the 
exponential type potential. 

Recently, Randall and Sundrum |l5|, [l6| have pointed out that a scenario with an infinite fifth dimension in the 
presence of a brane can generate a theory of gravity which mimics purely four-dimensional gravity, both with respect 
to the classical gravitational potential and with respect to gravitational radiation. The gravitational self-couplings are 
not significantly modified in this model. This result has been obtained from the study of a single 3-brane embedded 
in five dimensions, with the 5D metric given by ds 2 = e^^ v ^rj tlu dx fl dx u + dy 2 , which can produce a large hierarchy 
between the scale of particle physics and gravity due to the appearance of the warp factor. Even if the fifth dimension 
is uncompactified, standard 4D gravity is reproduced on the brane. In contrast to the compactified case, this follows 
because the near-brane geometry traps the massless graviton. Hence this model allows the presence of large or even 
infinite non-compact extra dimensions. Our brane is identified to a domain wall in a 5-dimcnsional anti-de Sitter 
space-time. 

The Randall-Sundrum (RS) model was inspired by superstring theory. The ten-dimensional Eg x Eg heterotic 
string theory, which contains the standard model of elementary particle, could be a promising candidate for the 
description of the real Universe. This theory is connected with an eleven-dimensional theory compactified on the 
orbifold R w x S 1 /Z 2 @. In this model we have two separated ten-dimensional manifolds. 

The static RS solution has been extended to time-dependent solutions and their cosmological properties have been 
extensively studied fL8|-[p7| (for a review of dynamics and geometry of brane universes see p8|]). 

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are 
confined on the 3-brane in a 5-dimensional space-time with ^-symmetry have been obtained, by using an elegant 
geometric approach, by Shiromizu, Maeda and Sasaki p9| , |30(| . The correct signature for gravity is provided by the 
brane with positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is 
required to be spatially homogeneous. The electric part of the 5-dimensional Weyl tensor E]j gives the leading order 
corrections to the conventional Einstei n eq uations on the brane. The effect of the dilaton field in the bulk can also 



be taken into account in this approach 1 3 1 



The linearized perturbation equations in the generalized RS model have been obtained, by using the covariant 
nonlinear dynamical equations for the gravitational and matter fields on the brane, by Maartens j32j. The behavior of 
an anisotropic Bianchi type I brane- world in the presence of inflationary scalar fields has been considered by Maartens, 
Sahni and Saini f33j. A systematic analysis, using dynamical systems techniques, of the qualitative behavior of the 
FRW, Bianchi type I and V cosmological models in the RS brane world scenario, with matter on the brane obeying a 
barotropic equation of state has been performed by Campos and Sopuerta |?4|, In particular, they constructed the 
state spaces for these models and discussed what new critical points appear, the occurrence of bifurcations and the 
dynamics of the anisotropy. The general exact solution of the field equations for an anisotropic brane with Bianchi 
type I and V geometry, with perfect fluid and scalar fields as matter sources has been found in |36|] . Expanding Bianchi 
type I and V brane-worlds always isotropize, although there could be intermediate stages in which the anisotropy 
grows. In spatially homogeneous brane world cosmological models the initial singularity is isotropic and hence the 
initial conditions problem is solved [ [37| . Consequently, these models do not exhibit Mixmaster or chaotic- like behavior 
close to the initial singularity j38|. 

Realistic brane-world cosmological models require the consideration of more general matter sources to describe the 
evolution and dynamics of the very early Universe. Hence the effect of the bulk viscosity of the matter on the brane 
have been analyzed in [ |39| . Limits on the initial anisotropy induced by the 5-dimensional Kaluza-Klein graviton 
stresses by using the CMB anisotropies have been obtained by Barrow and Maartens p0[. Anisotropic Bianchi type 



3 



I brane- worlds with a pure magnetic field and a perfect fluid have also been analyzed . 

It is the purpose of the present paper to investigate the effects of the rotational perturbations on a brane world with 
FRW type geometry. A similar case was also considered in . Assuming the rotation is slow and by keeping only the 
first order rotational terms in the field equations, a rotation equation describing the time and space evolution of the 
metric perturbations is obtained. This equation also contains the angular velocity of the matter rotating on the brane. 
By assuming that the metric perturbation and the matter angular velocity are separable functions of the variables t 
and r, the mathematical consistency of the rotation equation leads to some restrictions on the functional form of the 
angular velocity. In particular a class of solutions of the field equations leads to a barotropic brane- world cosmological 
model, with a non-linear pressure-energy density dependence. But generally, similar to the general relativistic case, 
the rotational perturbations will rapidly decay due to the expansion of the Universe. However, this general result is 
valid only in the presence of the dark energy term, describing the influence of the five-dimensional bulk on the brane. 
If this term is zero, rotational perturbations in a very high density (stiff) cosmological fluid decay only in the so-called 
case of the "perfect dragging" . 

The present paper is organized as follows. The field equations for a slowly rotating brane- world are written down 
and the basic rotation equation is obtained in Section II. The integrable cases of the rotation equation are considered 
in Section III. A brane world cosmological model, derived from the mathematical consistency requirement of the 
rotation equation, is obtained in Section IV. In Section V we discuss and conclude our results. 



II. GEOMETRY, BRANE- WORLD FIELD EQUATIONS AND CONSEQUENCES 

In the 5D space-time the brane-wor Id is located as Y(X I ) = 0, where X 1 , I = 0,1,2,3,4 are 5-dimensional 
coordinates. The effective action in five dimensions is |31(1 



S = J d 5 X^g~ 5 (j^R 5 - A 5 ) + jf ^ d*X^ ^Jf± - A + L-atter^ 



(2) 



with k\ = 8irG5 the 5-dimensional gravitational coupling constant and where x M , /i = 0,1,2,3 are the induced 4- 
dimensional brane world coordinates. We use a system of units so that the speed of light c = 1. i?5 is the 5D intrinsic 
curvature in the bulk and K is the extrinsic curvature on either side of the brane. 

On the 5-dimensional space-time (the bulk), with the negative vacuum energy A5 as only source of the gravitational 
field the Einstein field equations are given by 

Gu = klT u , Tu = -A 5giJ + S(Y) [-Xg u + Tfr attor ] , (3) 

In this space-time a brane is a fixed point of the Z2 symmetry. In the following capital Latin indices run in the range 
0, ...,4 while Greek indices take the values 0, 3. 

Assuming a metric of the form ds 2 = (?ijnj + gi^dx 1 dx J , with njdx 1 = d\ the unit normal to the x = constant 
hypersurfaces and gu the induced metric on \ = constant hypersurfaces, the effective four-dimensional gravitational 
equations on the brane take the form j2^, |30| : 

G>„ = —hg^v + k\T^ v + k^S^ — -E^jy, (4) 

where 

= -^TT^ - -^ a T va + ± 9iat (3T^ 3 T aP - T 2 ) , (5) 

and A = fc 2 (A 5 + fc 2 A 2 /6)/2, k\ = fcfA/6 and Eu = CiAjBn A n B . Ciajb is the 5-dimensional Weyl tensor in the bulk 
and A is the vacuum energy on the brane. T^ v is the matter energy-momentum tensor on the brane and T = ^ is 
the trace of the energy-momentum tensor. 

For any matter fields (scalar field, perfect fluids, kinetic gases, dissipative fluids etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as 

T^u = {p + p)u^u, y + ph^ + n^v + 2q( p u l/ ) . (6) 

The decomposition is irreducible for any chosen 4- velocity . Here p and p are the energy density and isotropic 
pressure, and h^ v = g^ v + u^u^ projects orthogonal to u^. The energy flux obeys = q< M >, and the anisotropic 
stress obeys 7r M „ = 7r </JV >, where angular brackets denote the projected, symmetric and tracefree part: 



V <ti> = VK, W <flu> 



W a p. (7) 
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The symmetric properties of E^ v imply that in general we can decompose it irreducibly with respect to a chosen 
4- velocity held it' 1 as 



E„ 



6 
Afcf 



U 



(8) 



where U is a scalar, a spatial vector and V^ u a spatial, symmetric and trace-free tensor. For a FRW model 
= V^u — p5[ and hence the only non-zero contribution from the 5-dimensional Weyl tensor from the bulk is 
given by the scalar term hi. 

The Einstein equation in the bulk imply the conservation of the energy momentum tensor of the matter on the 
brane, 



T^-v | x =o— 0. 

The rotationally perturbed metric can be expressed in terms of the usual coordinates in the form [|43[ 

dr 2 



ds z = -dt z + a 2 (t) 



1 — kr 2 



(<ffl 



2 ' sin 2 



2n(t,r)a 2 (t)r 2 sin 2 9 dtdip, 



(9) 



(10) 



where Q(t, r) is the metric rotation function. Although f2 plays a role in the "dragging" of local inertial frames, it is 
not the angular velocity of these frames, except for the special case when it coincides with the angular velocity of the 
matter fields, k = 1 corresponds to closed Universes, with < r < 1. k = — 1 corresponds to open Universes, while 
the case k — describes a flat geometry, where the range of r is < r < oo. In all models, the time- like variable t 
ranges from to oo. 

For the matter energy-momentum tensor on the brane we restrict our analysis to the case of the perfect fluid 
energy-momentum tensor, 



T" v = {p + p)u^u v +pg flv . 



(11) 



The components of the four- velocity vector are u — 1, u 1 



and u 3 = ui(t, r). uj — dip/dt is the angular velocity 



of the matter distribution. Consequently, for the rotating brane the energy-momentum tensor has a supplementary 
component 



T 03 = {[fi(t, r) - w(t, r)] p - w(t, r)p} r 2 a 2 (t) sin 2 



(12) 



We assume that rotation is sufficiently slow so that deviations from spherical symmetry can be neglected. Then to 
hrst order in Q the gravitational and held equations become 
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(13) 

(14) 

(15) 
(16) 



3«^M + ^M = 0, (17) 
a or atar 

(1 - kr 2 ) d2i l {t ; r) + (-- Mr) _ 4 [n(t, r) - w(t, r)] (k - ad + d 2 ) = 0. (18) 



dr 2 \r J dr 

The last two equations follows from the Ru and i?03 components of the held equations, respectively. Generally, we 
shall assume that the thermodynamic pressure p and the energy density p are related by a barotropic equation of 
state p = p(p). 

From a mathematical point of view the held equations ©-(U) and ©-fljj) represent a system of five equations 
in five unknowns a(t), p(t),U, f2(i, r) and u>(t, r) (the Bianchi identity (|15|) is a consequence of the held equations and 
the pressure can be eliminated via the equation of state of the matter on the brane). Therefore to hnd the general 
solution of this system is a well-posed problem. Since the Bianchi identity and the evolution equation for U can be 
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generally integrated, giving the energy density p and the dark matter term as some functions of the scale factor a, 
the basic field equations describing the first order rotational perturbations of brane world cosmologies are Eqs. jl^ ) 
and (|l7|)-(|l8|), three equations in three unknowns a(t),Q(t,r) and uj(t,r). 

Mathematically, Eqs. (0)-((l^) represent a system of partial differential equations. To obtain the solution of these 
equations we shall use the standard technique of the separation of variables [H. Therefore, we generally assume that 
both the rotation metric function and the angular velocity of the matter on the brane are the product of two functions, 
the first one depending on the cosmological time t only and the second one on the radial variable r only. When these 



representations of the unknown functions are substituted back to Eqs. (17)-([L8|), the field equations become identities 
in the independent variable, which can be separated in two independent equations, giving the time evolution and the 
spatial behavior of the considered physical parameters. If the values of the separation constants are properly chosen, 
then the combination of the solutions of the separated equations gives the solution of the initial partial differential 
equation for any boundary or initial conditions . 

Eq. ( |l6| ) can be immediately integrated to give the following general expression for the "dark energy" U: 

U = % (19) 

with Uq > a constant of integration. 

The temporal and spatial dependence of Q(t,r) is determined by Eqs. ( |l7| ) and (|l8|). The general solution of Eq. 
( |l7| ) can be immediately found and is given by 

where A{r) is a function to be determined from the field equations and an arbitrary time dependent function has 
been set to zero, without altering the physical structure of the model. Substituting Eq. (|2C]) into Eq. (|l^) gives the 
following equation describing the evolution of the rotational perturbations on the brane: 

(1 _ fcr 2 )^P +(-- 5kr) - 4 \A(r) - w(t, r)a 3 } (k -ad + a 2 ) = 0. (21) 

dr z \r ) dr 

By introducing a new variable r\ = r 2 , Eq. ( ^i"| ) can be transformed into the following form: 

v{1 k ^ )( ^r + {l~ 3kT] ) = [A{v) ~ v)a3] {k ~ a&+ h2) - (22) 

Eq. (p2|), the basic rotation equation, governs the evolution of the first order rotational perturbations of brane 
world models. For a given equation of state of the matter, a(t) is obtained from the field equations (|l3|)-(|l6|). Hence 
Eq. (||^) places some restrictions on the possible form of the angular velocity of the matter u(t, r). Indeed, since the 
left hand side of the Eq. (|2^) is a function of r alone, the right hand side must be either a function of r alone or a 
function of t alone (and hence set equal to a constant). As a first consequence of the mathematical structure of Eq. 
( ^2| ) it follows that the angular velocity of the rotating brane world must also be a separable function and we assume 
it to be of the form u)(t, r) — f(t)g(r), with f(t) and g{r) functions to be determined from Eq. (0). 



III. CONSISTENCY CONDITIONS FOR THE ROTATION EQUATION 

The first and simplest case in which Eq. (|22|) has a solution corresponds to the case u)(t, r) — Q(t, r). Physically, this 
situation corresponds to the so-called "perfect dragging", in which the function Q(t,r) appearing in the rotationally 
perturbed line element Eq. is the angular velocity of the rotating matter. Then the spatial evolution of the 

angular velocity is determined from the equation 

, (1 - M ^ + (|-3^^) =0 . ,23, 

Eq. (^3| ) has the following solutions: 



dr] 2 \2 J dr] 



C[ x \l-r 2 ) 1 ' 2 (2r- 1 +r- 3 ) k = +1, 

A(r) = <j C^ 1) (l + r 2 ) 1 /2( 2r -i_r- 3 ) +C^ 1) , k = -1, (24) 
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FIG. 1: Dynamics of the temporal part of the angular velocity f(a) (in a logarithmic scale) as a function of the scale factor a 
for different equations of state of the brane matter: 7 = 2 (solid curve), 7 = 4/3 (dotted curve) and 7 = 1 (dashed curve). We 
have normalized the parameters so that po = X, 8U0 = k%p% and 2C = fcfpo- 



(i) (i) 

where C{ and C\ , i = — 1, 0, 1 are arbitrary constants of integration. 

These solutions are not regular at the origin. They could be used in regions away from the origin and joined 
continuously to other solutions which are regular at the origin, which could be obtained by taking into account second 
order or higher terms in f2 in the field equations. The temporal behavior of the angular velocity is given by a~ 3 (t). 
Thus for an expanding Universe this type of dependence guarantees the decay of the rotational perturbations in the 
limit of large cosmological times. 

The second case in which the rotation equation can be solved corresponds to the choice g(r) = A(r). In this case 
both u>(t, r) and f2(i, r) have the same spatial dependence, but their time dependence is different. Therefore Eq. j2^ ) 
can be decoupled into the following two simple ordinary differential equations: 

*-")^ + (!- sfa ')^- Ci4W = 0, (25) 

[l-/(i)a 3 ] (fc-aa + a 2 ) = C, (26) 

where C is a separation constant. 

Consider first the time dependence of the angular velocity. With the use of the field equations Eqs. Eq. 
© gives 

« 2 [kl{p + P) + (fcfMV 2 + {kl/\)p P + (&VAfc 4 2 )q- 4 ] - 2C 
fW a 5 ikl(p + p) + (kj/\)p2 + (kl/\)pp+(8Uo/\kj)a- 4 } ' 1 ' 

In the very early stages of the evolution of the Universe it is natural to assume that the very high density matter 
obeys a barotropic equation of state of the form p = (7 — l)p, with 7 a constant and 1<7<2. 7 = 2 cor resp onds 
to the extreme limit of very high densities in which the speed of sound equals the speed of light (stiff case) |4q] . For 
a barotropic equation of state the conservation equation Eq. ( |l5| ) can immediately be integrated to give p = po a ~ 37 - 

Therefore for a cosmological fluid obeying a barotropic equation of state the time dependence of the angular velocity 
of the slowly rotating brane world is given by 

= IPokja 2 [a 3 i + (8ZV 7 Afc|p )a^- 4 + ( Po /\)} - 2Ca^ 
JW 7 Pofc 4 2 « 5 [« 37 + (8ZV7A£»a 6 ^ 4 + (Po/A)] ' [ ' 

The variation of the function f(a) for different values of 7 is represented in Fig. 1. 

For all values of 7 and for an expanding Universe the rotational perturbations tend, in the large time limit, to zero. 
This result is also independent on the numerical value of the separation constant C . 

In the particular case of the vanishing dark matter, Uq = 0, and in the limit of large times f(a) tends to the limit 



(29) 
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Therefore in the absence of dark energy the rotational perturbations on the brane would decay in time only if the 
barotropic fluid satisfies the condition 7 < 5/3. In particular this condition is not satisfied by the stiff cosmological 
fluid with 7 = 2. Hence in this case the observational evidence of a non-rotating Universe requires the condition 
C = and, consequently, a "perfect dragging" of the cosmological fluid. But in a radiation fluid with 7 = 4/3 the 
rotational perturbations will decay in time. 

For this model the spatial dependence of the angular velocity is described by Eq. (p5|). For k = +1 this equation is 
of the form x(l — x)F xx + [8 — (1 + a + (3)x\ F x — a(3F = 0, (the hypergeometric equation), with S 0, 1, 2... and a, (3 
are constants. The general solution of the hypergeometric equation is given by F — C±F(a, (3; <5; x) + C2X 1 ~ S F(1 — 8 + 
a, 1 — 5 + /?; 2 — S\ x), where C\ and C2 are arbitrary constants of integration, and F(a, (3; 8; x) — X^^Lo ^"nu}) x "' 
where (a)o = 1, (ct)n = r( a + n )/^( a ) — a i a + i)—(a + n — 1), n — 1,2, ... [Q. The radius of convergence of the 
solution is unity and if one of the constants a, (3, 8 — a, 8 — (3 is a negative integer, then the series terminates p6[ . 

For Eq. (|2^) 5 = 5/2 and the constants a and (3 must satisfy the conditions a + (3 = 2 and a/3 — C. The separation 
constant C must satisfy the condition C < 1. 

Therefore the general solution of Eq. (|||) is given by 



n— V2/n n— V ^/7l 



Since the second term is not regular at r = we must take = 0. Therefore the resulting solution is also defined 
and is regular at the origin of the radial coordinate. 

For open models k — —1 and the solution of Eq. fl25| ) can be obtained from the previous one by replacing rj by — 77 
and the separation constant C by — C. These solutions are convergent only for r\ < 1. For solutions regular outside 
rj = 1, one must consider solutions found in the neighborhood of the remaining two singular points of the differential 
equation, namely 1 and 00 |f46|. 



For k = Eq. (Elh gives 



with the general solution given by (wl 



d 2 A(r) 4 dA(r) 
dr 2 r dr 



AC'A(r) = 0, (31) 



g(r) = A(r) = 1 



(0)„-2VCr 



e 



A: = 0, (32) 



and and Cj '* arbitrary constants. Hence we have obtained in a closed form all the possible spatial distribution 



functions corresponding to the time dependence Eq. < \27\ ) of the angular velocity of the slowly rotating brane world. 

The third case in which the variable in Eq. ( pl| ) can be separated is given by assuming that ui(t, r) = G(r)/a 3 (t) 
and the time only dependent term in the right hand side of the equation is a constant. Therefore we obtain: 

- ^-w 1 + G " 3kv ) d ~^r~ K [m ~ G{v)] = °' (33) 

k-ad + h 2 = K, (34) 

where K is a separation constant. Generally Eq. ([33]) cannot be solved, due to a lack of exact knowledge of the 
mathematical form of the function G(r?). In the particular case G(r/) = Gq — constant, the solutions of Eq. ( |33| ) can 
be obtained by the substitution A(r) — > A(r) + Gq, where A(r) is any of the solutions of the homogeneous equation 
Eq.©. 

Finally, we consider the possibility of the existence of an uniformly rotating brane world, with lu — luq — constant. 
In this case the rotation equation takes the form 

r/(l - k v ) d2 f^ Skrj) " AMQt ~aa + a 2 ) = -^(k -aa + a 2 ). (35) 

dr\ \2 / dr] 



Therefore the mathematical consistency of Eq. (35) requires either A(rj) — constant, implying a — constant or 
(k — ad + a 2 ) — constant and a 3 (k — ad + a 2 ) = constant, also leading to a = constant. Hence uniform rotation is 
possible only for a static brane, the expansion of the Universe making the angular velocity a complicated function of 
the temporal and spatial coordinates. 
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IV. AN INFLATIONARY BRANE WORLD COSMOLOGICAL MODEL 



In the previous Section we have shown that the third condition of integrability of the rotation equation leads to the 
Eq. (|34|), describing the evolution of the scale factor of the slowly rotating brane world. Hence in this case the time 
evolution of the brane Universe is strongly correlated with its rotational properties. 

With the help of the substitutions a = u, u 2 = v, Eq. (^) can be transformed into a first order linear differential 
equation of the form 

dv 2 2 

= -v--B, (36) 



a a 



da 

where we denoted B — K — k. The general solution of Eq. j3^ ) is v = HqO 2 + B, with Hq a positive-definite arbitrary 
constant of integration. We will assume Hq > for later convenience. 
Therefore the general solution of Eq. ([34]) is given by: 

VBHq 1 sinh [H (t - t )} , B > 0, 

y/lBlH^cashiHoit-to)], B<0, (37) 
e H (t-t ) ) B = 0, 

where t is a constant of integration. 

In this class of models the Hubble parameter is given by 

H coth.[H (t-t )], B > 0, 
H = { H Q tanh [H (t - t )} , B < 0, (38) 
Hq = constant, B = 0. 



The energy density and pressure of the matter on the brane follow from the field equations Eqs. (|13J) and (|14[), with 
the energy-density of the cosmological fluid given by 




\k 2 A a 2 \ 2 kj a 4 + \kl \k\ j ' { ' 

In order to obtain Eq. ( j39|) we have used the identity 3d 2 /a 2 + 3k/a 2 = 3K/a 2 + 3Hq. 

Due to the presence of dark matter term, the energy density of the brane world has a maximum, dp/da = 0, 
corresponding to a value of the scale factor given by a max = (2/ki) ^JUq/KX. For this value of a the energy 
density satisfies the condition d 2 p/da 2 < 0. The maximum value of the energy density is given by /O max = 

A '3K 2 j AUq + 1 — 2A/Afc 2 + &H^/\k 2 — lj. If U = and the scale factor a is a monotonically increasing function 

of time, then the energy density of the brane world is a monotonically decreasing function for all times. 

The condition of the non-negativity of the energy density p > can be reformulated as a condition on the scale 
factor of the form (A — 3i?Q)Afc|a 4 — 3XKk 2 a 2 + 6Uq < 0, inequality that is satisfied for all a if 3Hq > A and 
3\K 2 k 2 + 8(3Hq — A)Uq < 0. Since these two conditions cannot be satisfied simultaneously, it follows that generally 



the energy density is non-negative only for a finite time interval. For Uq = 0, p > for a € (0, \/3K/ (A — 3Hq) 
The general pressure-energy density relation (the equation of state) is given by 

= m/kpa- 2 - (8^ /Afc|)a- 4 -p- (p 2 /X) 
1 + p/X 

The condition of the non-negativity of the pressure requires the condition 

2K 1 8W 1 / p\ 

be satisfied for all a. From Eq. ([l3|) we obtain 

f p\ 6K 1 12% 1 2(3ff 2 -A) 

p V + x) = ^-^kJ^- p + — kl — ' (42) 
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and therefore the condition of the non-negativity of the pressure can be reformulated as the following condition which 
must be satisfied by the energy density of the matter on the brane: 

AK 1 4U Q 1 | 2(3g 2 -A) ^ 
k^ oP" \k^ cl k ^ 

Hence with the use of Eqs. (|37]) we obtain the following exact analytic representations for the energy density and 
pressure: 



p{t) = A { Van sinh- 2 [H (t - to)] - b sirm" 4 [H (t - to)] +cq-1\, U Ml. ! 44 ! 

k sinh" 2 [H (t - t Q )] - u sinh" 4 [H Q (t - t )] 



P(t) 



a sinh 2 [H (t - t )] - b sinh 4 [H (t - to)] + c Q 



\iJaosmh- 2 [H (t-to)]~b smh~ 4 [Ho(t-t Q )]+c Q -l\, B > 0, (45) 



p(t) = X^ao cosh" 2 [H (t - t )] - b cosh" 4 [H (t - to)] ~+cq - 1 J , B < 0, (46) 

fc cosh" 2 [JJ (* - t a )] - uo cosh- 4 [H (t - i )] 
= j = 
y a cosh 2 [if (t - to)] - b cosh 4 [ff (t - f )] + c 

a cosh- 2 [J? (i -to)] - fcocosh" 4 [H (t-t )] + c - ll , B < 0, (47) 



K*) = A |y^ ex P[- 2 ^o(t-to)]-^e X p[-4ffo(t-*o)]+co-l|, B = 0, (48) 

= (2^/fc 2 ) exp [-2ff (t - t )] - (8^o/Afc 2 ) exp [-4ffo(t - t )] 

V(6^/Afc 2 ) exp [-2H (t - to)] - (12^o/A 2 fc|) cxp [-AH Q {t - to)] + c 

- A | ^{&K/\kD exp [-2H {t - to)] - (l2U /\ 2 kj) exp [-4H (t - to)] + c - l| , 5 = 0, (49) 

where we denoted a = GKHg/Xh^B], 6 = 12W ^/A 2 ^|S| 2 , c = 1 - 2A/Afc| + 6#£/Afc|, k = 2KH^/kl\B\ and 
mo = 8W ^/Afcl|B| 2 . 

An important observational quantity is the deceleration parameter, defined as q = (d/ dt)H~ l — 1, given by 

-tanh 2 [H (t - t )] , £ > 0, 

-coth 2 [ff (t-to)], B<0, (50) 
-1, B = 0. 

There are two distinct types of behavior for the brane-world type cosmological models described by Eqs. (f44|)- 
([49|). For the first model, with the scale factor having a singularity at the origin, with an appropriate choice of the 
parameters, the energy density can be scaled to zero at the beginning of the cosmological evolution, corresponding to 
t = to. Then, for times t > to, the energy density is an increasing function of time and reaches a maximum value at 
t = t max . For time intervals t > t max the energy density is a decreasing function of time. Hence this model can be 
used to describe matter creation on the brane, a phenomenon which is entirely due to the presence of the dark matter 
term. If the dark energy term IA — 0, then, due to the singular behavior of the scale factor, the initial energy density 
and pressure of the cosmological fluid are infinite at the beginning of the cosmological evolution, p,p — > oo for a — » 0. 

For the second and third class of models, the scale factor, energy density and pressure are all finite at t = to- Hence 
in this case the brane is filled with an initial cosmological fluid. Depending on the numerical choice of the parameters, 
there are also two types of cosmological behavior, with the energy density reaching its maximum value at the initial 
moment or at a time t — t max > to- In the first case for times larger than the initial time, the energy density and the 
pressure are monotonically decreasing functions for all times t > to, while in the second case p and p are monotonically 
decreasing functions of time only for t > t max . 
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FIG. 2: Time evolution of the energy density p(t) for the brane world cosmological model, for the three classes of solu- 
tions: first class (o ~ sinh [Ho(t — to)]) (solid curve), second class (a ~ cosh [Ho(t — to)]) (dotted curve) and third class 
(a — exp [Ho(t — to)]) (dashed curve). We have normalized the parameters so that IK = k\, QK/\k\ — 1, 12Wo/A 2 fc| = 1 and 
2(3Hq — A) = Afc|. For the sake of presentation the curves have been rescaled with different factors. 
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FIG. 3: Time evolution of the pressure p(t) of the cosmological fluid for the brane world cosmological model, for the three 
classes of solutions: first class (a ~ sinh [Ho(t — to)]) (solid curve), second class (a ~ cosh [Ho(t — to)]) (dotted curve) and third 
class (a = exp [Ho(t — to)]) (dashed curve). We have normalized the parameters so that 2K = k\, §Kj\k\ — 1, 12Wo/A 2 fc| = 1 
and 2(3-ffo — A) = Afe|. For the sake of presentation the curves have been rescaled with different factors. 

The behavior of the energy density and pressure is represented, for all three cosmological models and for a particular 
choice of the numerical values of the parameters, in Figs. 2 and 3. 
The equation of state of matter is presented in Fig. 4. 

For the given choice of parameters, the pressure-energy density dependence is generally non-linear, with a depen- 
dence which can be approximated for small densities by a linear function, p ~ p. In this region the effects of the 
contribution from the extra-dimensions can be neglected. For high densities the non-linear effects become important, 
showing that the quadratic effects due to the effects of the five-dimensional bulk also modify the equation of state of 
the matter. 

An important condition that must be satisfied by any realistic equation of state of dense matter is the requirement 
that the speed of sound c s = {dp/dp) 1 / 2 be smaller or equal to the speed of light, c s < 1. The time variation of the 
speed of sound in the cosmological fluid with equation of state given by Eq. (BO) is represented in Fig. 5. 

For the given range of parameters and in the considered time interval, the condition c s < 1 is satisfied. During this 
period the speed of sound is a rapidly increasing function of time for the first solution and a slowly increasing function 
of time (almost a constant) for the second and third class of solutions. But for other choices of the numerical values 
of the constants or different time intervals this condition could be violated. 
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FIG. 4: Equation of state of the cosmological fluid for the brane world cosmological model, for the three classes of solu- 
tions: first class (a ~ sinh [Ho(t — to)]) (solid curve), second class (a ~ cosh [Ho(t — to)]) (dotted curve) and third class 
(a = exp [Ho(t — to)}) (dashed curve). We have normalized the parameters so that 2K — fcf , 6K/\k% = 1, 12Wo/A 2 fc| = 1 and 
2(3ff ( 2 — A) = Afcf . For the sake of presentation the curves have been rescaled with different factors. 
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FIG. 5: Time variation of the speed of sound c s = (dp/ dp) 1 / 2 in the cosmological fluid for the three classes of solutions: first class 
(a ~ sinh [Ho(t — to)]) (solid curve), second class (a ~ cosh [Ho(t — to)]) (dotted curve) and third class (a = exp [Ho(t — to)]) 
(dashed curve). We have normalized the parameters so that 2K = fc|, QK/\k\ = 1, 12Wo/A 2 A;| = 1 and 2(3i/ ( 2 — A) = \k\. 

Since q < for all times, the evolution of the models is purely inflationary. However, due to the corrections from 
the extra-dimensions, the initial period of the inflationary phase can be associated to an increase in the energy-density 
of the Universe. 

In all these cases, due to the rapid expansion of the brane-world, there is a fast time decay of the angular velocity 
uu ~ a~ 3 . Since in the large time limit the scale factor is of the form exp(H t), in this brane world model one obtains 
an exponential decay of the cosmic vorticity. 

V. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the evolution of the rotational perturbations in the framework of brane 
world cosmology. As a first step we have obtained the basic dynamical equation governing the temporal variation and 
spatial distribution of the metric perturbation function Q(t, r). This equation also contains the angular velocity of 
the slowly rotating matter. From the consistency condition of the rotation equation one can obtain some restrictions 
on the admissible mathematical form of u. As a general result we have shown that in the presence of dark energy, 
describing the effects of the bulk on the brane, rotational perturbations always decay for the slowly rotating brane. 
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But for a vanishing dark energy term, rotational perturbations do not decay for very high density stiff matter with 
7 = 2, except for the case of the perfect dragging. Hence the dark energy term could have play an important role in 
suppressing the vorticity of the very early Universe. 

In order to find the time evolution of the four-velocities of the particles on the brane we consider the dynamics 
of a test particle in the perturbed metric (|l0|). The equations of motion are du^/ds + T^ v \u u u x = 0, where are 
the components of the four-velocity and the Christoffcl symbols T^ v \ are computed from the metric. To simplify the 
calculations we consider only the first order corrections to the metric in f2 and assume that test particles have small 
velocities, thus retaining only terms which are linear in velocity. Consequently we obtain 

n du 1 2d i du 2 du 3 2d • 2d o . . 

u° = i, -j- = — u 1 , -jr = o, — = — n + n-— u 3 . (51) 



dt a dt dt a a 



By integrating Eqs. (|5l]) we find 



»' = j|r „* = „ s , »» = 4 J+ ^, (52, 

with Uq, i — 1,2,3 constants of integration. Uq can be set to zero without altering the physical structure. In the 
limit of large t, a(t) oo and we have u 1 = and u 3 = 0. At the time t — the test particle will have shifted in 
azimuth by Aip — J*°° u 3 dt, having an angular velocity u = ^j. Present-day observations impose a strong restriction 
on the numerical value of the angular velocity of the Universe of the form w(ioo) < 10 _15 years _1 . Hence in realistic 
cosmological models the angular velocity must tend to zero in the large time limit. 

The condition of the mathematical consistency of the rotation equation also leads to a specific brane world cosmo- 
logical model, with the equation of state of the matter given in a parametric form. In the high density regime the 
pressure-energy density dependence is nonlinear. The presence of the dark energy term has major implications on the 
time evolution of physical parameters. Due to the presence of U, the energy density and pressure of the matter on 
the brane have a maximum. The increase in the energy density from zero to a maximum value can model the energy 
transfer from the bulk to the zero-density brane, thus leading to the possibility of a phenomenological description of 
matter creation on the brane. If the dark energy term is zero, the singularity in the scale factor is associated to a 

singular behavior of p and p. In the large time limit, a — > oo and p — > A — 2A/Afc| + QH^/Xkl — lj , with the 

pressure p — > —p. Since for ordinary matter p cannot be smaller than zero, it follows that in this limit p must be 
zero, condition which leads to H = \/A/3 and a = exp(ifot)- Therefore the de Sitter solution is an attractor for this 
brane world model. Independently on the initial state, the Universe ends in an inflationary phase. 

On the other hand one must point out that the curvature of the brane, described by the parameter k = —1,0, +1 
does not play any significant role in the time evolution of this brane- world model. Since k is absorbed in the arbitrary 
constant B (which in fact is the most important parameter of the model), the dynamics of the matter on the brane is 
independent on the three-dimensional geometry of the Universe. However, the knowledge of the spatial distribution 
of the angular velocity could, at least in principle, allow for the determination of the separation constant K . If 
the numerical value and sign of K would be known, then the effects of the curvature of the brane on the temporal 
dynamics could also be obtained. But due to the exponential increase in the scale factor, the Universe ends in a 
non-rotating state with a flat geometry. 

Since the equation of state of matter is unusual (even in the extreme limit of high densities the equation of state of 
ordinary matter is still linear with p = p), it is natural to assume that in this cosmological model the initial matter 
content of the Universe on the brane consists from a field rather than ordinary matter. The best candidate is a scalar 

5 2 12 

field <f>, with p = p0 = ^- + V(4>) > and p — p^, = \ — ^('Wj where V(<p) > is the self-interaction potential. Then 
the time evolution of the scalar field can be immediately obtained from <j)(t) = \/2(p + p), while the potential is given 
by V(t) = (p — p)/2. With the use of Eqs.([l4|)-([49|) one can easily obtain the time dependence of the scalar field and 
scalar field potential for each class of solutions. The scalar field drives the Universe into an inflationary era. In this 
case solutions with negative pressure are also allowed. In the limit of large times, since the field satisfies the equation 
of state p + p = 0, we obtain = and <j> — constant. For the scalar field potential we find V — * p, and, in the de 
Sitter limit, it follows V = 0. Hence V does not give any contribution to the cosmological constant A. 
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